The dominant imprint of Rossby waves in the climate network 
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The connectivity pattern of networks, which are based on a correlation between ground level 
temperature time series, shows a dominant dense stripe of links in the southern ocean. We show 
that statistical categorization of these links yields a clear association with the pattern of an atmo- 
spheric Rossby wave, one of the major mechanisms associated with the weather system and with 
planetary scale energy transport. It is shown that alternating densities of negative and positive links 
(correlations) are arranged in half Rossby wave distances around 3,500 km, 7,000 km and 10,000 
km and are aligned with the expected direction of energy flow, distribution of time delays and the 
seasonality of these waves. It is also shown that long distance links (i.e., of distances larger than 
2,000 km) that are associated with Rossby waves are the most dominant in the climate network. 
Climate networks may thus be used as an efficient new way to detect and analyze Rossby waves, 
based on reliable and available ground level measurements, in addition to the frequently used 300 
hPa reanalysis meridional wind data. 
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Networks have become an important tool for analyzing 
technological and natural systems jl|-|3 • Examples range 
from social relations [5[, biochemical interactions [6|-[8|, 
information flow through the world wide web [9|, physi- 
ological activities ^3 5 ^^^ ^^^ mitigation of attacks on 
transportation infrastructures Hi]. It was suggested in 
past years that climate variables, like temperature and 
geopotential height, can be viewed as a climate net- 
work [12-19] . In this representation, different regions 
of the world are regarded as nodes of the network, and 
links of the network represent communications between 
different locations via, e.g., heat and material exchange. 
A multitude of statistical analysis methods are often used 
to capture major variability patterns in climate time se- 
ries |2Q|-[22| , where the correlation matrix plays an impor- 
tant part. The climate network provides a complemen- 
tary tool to study the statistical properties of the climate 
system. 

The climate network often has very stron g li nks which 
are caused by a proximity (distance) effect pjj . Namely, 
pairs of sites close to each other below some threshold dis- 
tance (typically 2,000 km) are often strongly positively 
correlated. A significant fraction of the stable network 
structure may be associated with the proximity effect. 
It is hence common to analyze climate time series also 
based on negative correlations (e.g., |23l-[25j), which usu- 
ally represent more interesting remote interactions, called 
teleconnections. 

However, recent studies have not distinguished be- 
tween positive and negative correlations in climate net- 
works [iqI, fl2l - [l9| . In particular, it is apparent, based on 
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previous studies, that a large fraction of the links in the 
climate network resides in the southern ocean |13Lll4JJl7|. 
This fraction may include (beyond the proximity effect 
distance) both negative and positive correlations. 

Here we analyze separately the negative and positive 
correlations of the climate network. We show that these 
links alternate, as a function of distance, between neg- 
ative and positive, consistent with a wave pattern. We 
find that the time delay associated with these links in- 
creases from one to five days as a function of the distance 
between the nodes (the length of the links). We also ana- 
lyze the typical length scale of the links, their seasonality, 
and the geographical structure of the climate network; 
all of these are found here to be consistent with atmo- 
spheric Rossby waves, one of the most efficient climate 
mechanisms of planetary-scale energy transfer [26| . Stud- 
ies of atmospheric Rossby waves are usually based on the 
300 hPa meridional wind velocity reanalysis data |27l-[29j . 
Here we show that it is possible to uncover the character- 
istics of Rossby waves using more common and reliable 
surface data, like surface air temperature. We find that 
Rossby waves dominate the climate network, an observa- 
tion that, surprisingly, has not been previously reported. 

Here, we analyze the daily data of air temperature, 
sea level pressure, geopotential height, and meridional 
velocity fields. Specifically, we analyze a network of 
726 nodes around the globe (see [13 or small dots in 
Fig. 4) from the NCEP/NCAR reanalysis I grid ^,^. 
Below, we mainly focus on the surface temperature 
field as one of the most common and reliable types 
of data. For each node (i.e., longitude/latitude grid 
point) of the network, daily values within the period 
1948-2010 are used, from which we extract anomaly val- 
ues. Specifically, if we take the record of a given site 
in the grid to be T^((i), where y is the year and d is 
the day (from 1 to 365), then the filtered record is de- 
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FIG. 1: (Color online) (a) An example of a seasonally de- 
trended surface temperature daily time series from 45° S, 
52.5°E for the year 1948. (b) Same as (a) for 45°S, 97.5°E. (c) 
The cross-correlation function between the time series shown 
in (a) and (b). (d) Same as (c) but when the time series of 
the two sites are of different years ("shuffling" method). 



noted by T^(d) = fy{d) - j^T^y Ty{d), where TV is the 
number of years available in the record. We also de- 

fine e,(d) = [Ts{d) - {Ts{d))]/(^{Ts{d) - {Ts{d))f)"\ 

where (• • • ) is the average of the time series. 

The link between each pair of sites on the grid, si 
and 52, is calculated as the cross-correlation function 
X|^ ,^(r > 0) = (e|^(^)e|^(d + r)), where r is the time 
lag and X|^ ^^{t) = XJ^ si(~'^)- We define the time lag, 
r*, at which XJ^ s^^^) ^^ maximal (or minimal), as the 
time delay of a pair si, S2. When si is to the west of S2 
and the time lag is positive, the link direction is to the 
east. We distinguish between positive and negative link 
weights as follows 
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MAX(X|^,,J-MEAN(X|'^,,J 
STD(X|'„.,) 



STD(X|'„,,) 
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(2) 



where MAX and MIN are the maximum and minimum 
values of the cross-correlation function, MEAN and STD 
are the mean and standard deviation, and the super- 
script y denotes a specific year. Typical time series and 
their cross-correlation functions are shown in Fig. [TJ In 
this sample, the absolute value of a minimal (negative) 
cross-correlation function is much larger than the maxi- 
mal value. 

We implemented the above procedure for winter and 
summer time series, separately. Specifically, we calculate 
the cross-correlation functions of each pair of sites, where 
the data ranges from May 1st to Aug. 31st (123 days) for 
southern hemisphere (SH) winter (northern hemisphere. 



(NH) summer) or from Nov. 1st to Feb. 28th (120 days) 
for SH summer (NH winter). We choose a maximal time 
lag of 72 days and a time period of ^ G [1948, 2010]. The 
correlation coefficient and the link weight are based on 12 
months of data, in order to have sufficient statistics. This 
is done by "gluing" together three consecutive winters 
(summers), such that the total number of months is 12. 
We analyze a near surface (1000 hPa) temperature 
time series. First, we focus on the properties of link 
weight around the globe (Fig. 2) for the months Nov. 
to Feb.. To identify the significant links we apply a shuf- 
fling procedure in which the order of the years is shuf- 
fled while the order within each year remains unchanged. 
This shuffling scheme is aimed at preserving all the sta- 
tistical quantities of the data, such as the distribution of 
values, and their autocorrelation properties, but omitting 
the physical dependence between different nodes. Fig. 
2(a) and (b) depict the link weight statistics for the real 
and shuffled data. High negative mean link weight values 
exist in the probability density function (PDF) of the real 
data but are missing in the shuffled data, and therefore 
are not likely to occur by chance. Moreover, high vari- 
ability (std.) during different years of the time delays 
of links r* is also a signature of random behavior [10| 
(Fig. 2(b)). The differences between the distribution 
of real data and shuffled data indicate that many sig- 
nificant negative links exist in the climate network (Fig. 
2(a),(b)). We obtain similar results for positive links and 
other fields (not shown). 

Next, we divide the world into three geographical 
zones, the SH (from 22.5°S to the south pole), the NH 
(from 22.5°N to the north pole) and the equator (between 
22.5°S and 22.5°N). We then calculate the (geographical 
and temporal) mean link weight {Ws-^^s^) ^s a function of 
the geographical distance of links, d. It is clear that in 
the SH, there is a preferred distance of ~ 3,500 km and 
a much weaker one of ~ 10,000 km (Fig. 2(c)). In the 
NH region, we find a similar, weaker dependence, while 
in the equatorial region, there is no preferred distance 
(Fig. 2(c)). These preferred distances may be associ- 
ated with atmospheric Rossby waves [26l-[29| , which have 
a wavelength of ~ 7,000 km and which are known to 
be pronounced in the SH, weaker in the NH and absent 
in the equatorial region [28|, [29| . The negative peaks at 
3,500 km and 10,000 km represent a 1/2 wavelength and 
a 3/2 wavelength of the observed Rossby wavelength. 

To further consolidate the association of the observed 
pattern in the climate network with Rossby waves, we 
compare the seasonality of this pattern with the known 
seasonal characteristics of Rossby waves. In Fig. 3(a)- 
(d), we plot the negative and positive weights of all 
SH links, for the winter and summer months separately. 
Each point represents an average link weight {W) over 
years versus its distance d. The negative weights, de- 
fined in Eq. ([2]), have a pronounced enhanced distri- 
bution of large weights for d ^ 3,500 km during both 
summer and winter months, while for the SH summer 
months (Nov. to Feb.), there is an additional preferred 
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FIG. 2: (Color online) (a) The probability density function (PDF) of the mean weight of negative links in the globe, in both 
real (blue solid line) and shuffled (red dashed line) data, (b) The dependence of the standard deviation of time delay r* on the 
mean weight (W) for all possible links, for real (blue) and shuffled (red) data, (c) The mean negative link weight as a function 
of distance d, for the SH (blue circles), NH (red squares) and the equator (pink diamonds) regions. The figure is based on the 
globe temperature records at 1000 hPa isobar for Nov. to Feb.. 
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FIG. 3: (Color online) The dependence of the weight of neg- 
ative links on the distance d in the SH during (a) winter and 
(b) summer months. (c),(d) Same as (a),(b) for positive links. 
The PDF of time delays r* of near negative (blue squares, d ^ 
3,500 km), far positive (red circles, d ^ 7,000 km), and far 
negative (pink diamonds, d ^ 10,000 km) weighted links, in 
the SH during (e) winter and (f) summer. The figure is based 
on the temperature field at 1000 hPa isobar. 



distance of ~ 10,000 km (Fig. 3 (b)), both in accor- 
dance with the 1/2 and 3/2 wavelengths of atmospheric 
Rossby waves. Around the full wavelength distance (i.e., 
d ^ 1^ 000 km), we find, as expected, an enhanced distri- 
bution of large positive weights (defined in Eq. ([1])) (Fig. 
3 (c),(d)). However, a larger abundance of strong waves 
(represented by links at 1/2, 1, and 3/2 wavelengths) dur- 



ing summer, in comparison to winter months, is clearly 
seen in Fig. 3 (a)-(d), in agreement with the clearer pat- 
tern of Rossby waves found during the SH summer in 
the SH [i^-il. The situation is similar in the NH (not 
shown), for the NH winter. In the distribution of pos- 
itive weights (Fig. 3(c),(d)), one clearly sees the links 
that emerge around d < 2,000 km due to the proximity 
effect mentioned in the introduction. 

Atmospheric Rossby waves have a characteristic group 
velocity, and we now estimate it based on the climate net- 
work results shown above. To achieve this, we divide d by 
r* for each link, under the assumption that r* is a good 
estimate for the underlying dynamical delay between the 
two sites (nodes) [32|- Since r* is only meaningful for 
links with weights above the background noise level (see 
Fig. 2(b), and ref. fig), we limit the analysis to links 
with weights \W\ > 2.8 [33|. Furthermore, to avoid links 
that are prone to the proximity effect (as such links can- 
not easily fit in the current Rossby wave interpretation) , 
when considering positively correlated links, we consider 
links with d > 5,000 km. Also, we constrain the near 
negative links to be within d G [2, 000/cm, 5, OOO/cm] and 
far negative links to have d > 8, 000 km. The PDFs of 
time delay r* of the near negative links, the far positive 
links and the far negative links are shown in Fig. 3(e),(f). 
By our convention, positive r* means an eastward energy 
flow, which is the typical case for most observed links. 
The results point to a time lag of r ^ 1 day for the near 
negative links (i.e., links in the first peak in Fig. 3(a)), 
r ^ 2-3 days for the far (beyond the proximity effect) 
positive links, and r ^ 3-4 days for the far negative links 
(i.e., belonging to the second peak in the scatter plot of 
Fig. 3(a,b)). Based on these numbers, our estimated 
group velocities are in the range [20 — 35]m/s, consistent 
with the range [23 — 32]m/s reported in previous stud- 
ies Hill. 

The climate network has a unique geographical struc- 
ture that can be compared with the geographical struc- 
ture of Rossby waves. Since the network is directed — 
positive r indicates eastward flow while negative r indi- 
cates westward flow — we distinguish between a link that 
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FIG. 4: (Color online) The out- (a) and in-degree (b) of the climate network structure from Nov. to Feb. (SH summer). 
The out- (c) and in-degree (d) of the climate network structure from May to Aug. (SH winter). The figure is based on the 
temperature record at 1000 hPa isobar. 



is pointing toward a node (where the number of links 
pointing to a specific node is referred to below as "in- 
degree"), or away from the node (referred to below as 
"out-degree") (see [13 )• Fig. 4 depicts the mean in- 
and out-degrees of each node, excluding the equatorial 
region that conforms with a pattern that is not related 
to our current discussion. The observed structure is con- 
sistent with the structure of Rossby waves [29]. First, the 
wave band in the SH from May to August (SH winter. 
Fig. 4(c), (d)) is broader than that of the SH summer 
(Nov. to Feb., Fig. 4 (a,b)). Second, the atmospheric 
Rossby wave structure in the NH summer is less pro- 
nounced. Third, the wave structure in the SH summer 
(Nov. to Feb.) lies on a band centered near 50° S. All the 
above characteristics are consistent with the properties 
of Rossby waves. 

Previous studies found the 300 hPa meridional veloc- 
ity field to be the most suitable for studying the charac- 
teristics of Rossby waves |27l-[29| . Our method captures 
the wave properties also using other fields at various al- 
titudes. In particular, we showed above that the wave 
pattern is clearly seen at a ground level (1000 hPa) tem- 
perature field, a more common and reliable variable. In 
Fig. 5, we compare, using the climate network technique, 
between the mean weight distributions of the negative 
links of the meridional velocity and temperature fields. 
We find that the pattern of the two fields is similar, al- 
though the meridional velocity yields larger weights. In 
addition, the meridional velocity yields a clearer pattern 
at the high altitude of 300 hPa, while the temperature 
field yields a clearer pattern at the ground level of 1000 
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FIG. 5: (Color online) The dependence of the weight of neg- 
ative links on the distance d in the SH during SH summer 
months for temperature (left panels) and meridional veloc- 
ity (right panels) for 300 hPa (upper panels) and 1,000 hPa 
(lower panels). 



hPa. 

In summary, we analyze the properties of the climate 
network by considering, separately, positive and negative 
correlations (links). The most dominant links in the cli- 
mate network with a geographical distance larger than 
2,000 km are found for distances of ~ 3,500 km, ~ 7,000 
km and ~ 10,000 km. These distances coincide with 
the 1/2,1 and 3/2 wavelengths of common atmospheric 



Rossby waves. Moreover, the time delays associated with 
these distances are in agreement with the direction of 
the energy flow and with the group velocity of the atmo- 
spheric Rossby waves. The pronounced length scales of 
the climate network, the dominance in the SH in com- 
parison with the NH, and the dominance during the SH 
summer in the SH are all consistent with the properties 
of atmospheric Rossby waves. All of these factors thus 
provide strong support for the association of the majority 
of the climate network far links with Rossby waves. 



ACKNOWLEDGMENTS 



The authors would like to acknowledge the support of 
the LING project (no. 289447) funded by the EC's Marie- 
Curie ITN program (FP7-PEOPLE-2011-ITN) and the 
Israel Science Foundation for financial support. We 
thank Edmund Chang and Nili Harnik for their assistance 
in understanding the characteristics of Rossby waves. 



[1] R. Albert, A-L. Barabasi, Rev. Mod. Phys., 74, 47 (2002). [19 

[2] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez and D-U. 

Hwang, Physics Reports, 424, 175 (2006). [20 

[3] M. E. J. Newman, Networks: An Introduction. Oxford 

Press, (2010). [21 

[4] R. Cohen and S. Havlin, Complex Networks: Structure, 

Robustness and Function, Cambridge University Press, [22 

(2010). 
[5] M. Girvan and M. E. J. Newman, Proc. Natl. Acad. Sci., 

U.S.A. 99, 7821 (2002). [23; 

[6] R. Guimera and L. A. N. Amaral, Nature, 433, 895 

(2005). [24 

[7] G. Palla, I. Derenyi, I. Farkas, and T. Vicsek, Nature, 

435, 814 (2005). [25 

[8] H. Jeong, B. Tombor, R. Albert, Z. N. Oltvai, A.-L. 

Barabasi, Nature, 407, 651 (2000). [26 

[9] B. A. Huberman, The Laws of the Web, MIT Press, Cam- [27 

bridge, MA (2001). [28^ 

[10] A. Bashan et al.. Nature Communications, 3, 702 (2012). 

[11] C. M. Schneider, A. A. Moreira, J. S. Andrade Jr, S. [29 

Havlin, H. J. Herrmann, Proc. Nati. Acad. Sci., U.S.A. [30 

108, 3838 (2011). 

[12] A. A. Tsonis, K. L. Swanson, and P. J. Roebber, Bull. [31 

Am. Meteorol. Soc, 87, 585 (2006). 
[13] K. Yamasaki, A. Gozolchiani, and S. Havlin, Phys. Rev. 

Lett., 100, 228501 (2008). 
[14] A. A. Tsonis and K. L. Swanson, Phys. Rev. Lett., 100, [32 

228502, (2008). 
[15] J. F. Donges, Y. Zou, N. Marwan and J. Kurths, Euro- [33 

phys. Lett., 87 48007, (2009). 
[16] Y. Berezin et al.. Nature Scientific Reports, 2, 666 (2012). [34 

[17] A. Gozolchiani, S. Havlin and K. Yamasaki, Phys. Rev. 

Lett., 107, 148501, (2011). 
[18] O. Guez et al., Eurphys. Lett., 98, 38006 (2012). 



J. F. Donges, Y. Zou, N. Marwan, and J. Kurths, Eur. 
Phys. J. Special Topics, 174, 157 (2009). 
H. V. Storch and F. W. Zwiers, Statistical Analysis in 
Climate Research, Cambridge University Press, (2003). 
L. F. Robinson, V. H. de la Pena and Y. Kushnir, The- 
oretical and Applied Climatology, 94, 215 (2008). 
P. Alpert and T. Sholokhman, Factor Separation in the 
Atmosphere: Applications and Future Prospects, Cam- 
bridge University Press, (2011). 

J. M. Wallace and D. S. Gutzler, Mon. Wea. Rev., 109, 
784 (1981). 

B. J. Hoskins and D. Karoly, J. Atmos. Sci., 38, 1179 
(1981). 

G. T. Walker, Quarterly Journal of the Royal Meteoro- 
logical Society 51, 337 (1925). 
E. K. M. Chang, J. Atmos. Sci., 62, 2467 (2005). 
E. K. M. Chang, J. Atmos. Sci., 50, 2038 (1993). 
E. K. M. Chang, and B. Y. Daniel, J. Atmos. Sci., 56, 
1708 (1999). 

E. K. M. Chang, J. Atmos. Sci., 56, 1729 (1999). 
Kalnay et al.. Bull. Amer. Meteor. S oc, 77, 437-47 0, 
(1996). http://www.esrl.noaa.gov/psd/data7iridded7] 
The original resolution of the reanalysis data is 2.5°. The 
number of nodes that were analyzed was normalized by 
the cosine of latitude, to keep the density of nodes equally 
distributed. 

L. Cimponeriu, M. Rosenblum and A. Pikovsky, Phys. 
Rev. E, 70, 046213 (2004) 

We have checked thresholds in the range 2.7 — 3 without 
significant qualitative changes in the results. 
E. H. Berbery and C. S. Vera, J. Atmos. Sci., 53, 468 
(1996). 



